Decay of the pseudoscalar glueball into scalar and pseudoscalar mesons 
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We study a chiral Lagrangian which describes the two- and three-body decays of a pseudoscalar 
glueball into scalar and pseudoscalar mesons. The various branching ratios are a parameter-free 
prediction of our approach. We compute the decay channels for a pseudoscalar glueball with a mass 
of 2.6 GeV, as predicted by Lattice QCD, which is in the reach of the PANDA experiment at the 
upcoming FAIR facility. For completeness, we also repeat the calculation for a glueball mass of 2.37 
GeV which corresponds to the mass of the resonance X(2370) measured in the BESIII experiment. 
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I. INTRODUCTION 



The fundamental symmetry underlying Quantum Chromodynamics (QCD), the theory of strong interactions, is 
the exact local SU(3) C color symmetry. As a consequence of the non-abelian nature of this symmetry the gauge fields 
of QCD, the gluons, are colored objects and therefore interact strongly with each other. Because of confinement, one 
expects that gluons can also form colorless, or 'white', states which are called glueballs. 

The first calculations of glueball masses were based on the bag- model approach [1| . Later on, the rapid improvement 
Q_[ of lattice QCD allowed for precise simulations of Yang-Mills theory, leading to a determination of the full glueball 
spectrum [2|]- However, in full QCD (i.e., gluons plus quarks) the mixing of glueball and quark-antiquark configurations 
with the same quantum number occurs, rendering the identification of the resonances listed in the Particle Data Group 
(PDG) more difficult. The search for states which are (predominantly) glueball represents an active experimental 
t-H ■ and theoretical area of research, see Ref. [H and refs. therein. The reason for these efforts is that a better understanding 
of the glueball properties would represent an important step in the comprehension of the non-perturbative behavior 
of QCD. However, although up to now some glueball candidates exist (see below), no state which is (predominantly) 
glueball has been unambiguously identified. 
\q ■ In general, a glueball state should fulfill two properties regarding its decays: it exhibits 'flavor blindness', because the 
gluons couple with the same strength to all quark flavors, and it is narrow, because QCD in the large- iV c limit shows 
that all glueball decay widths scale as N~ 2 , which should be compared to the N^ 1 scaling law for a quark-antiquark 
state. The lightest glueball state predicted by lattice QCD simulations is a scalar-isoscalar state (J PC = ++ ) with a 
mass of about 1.6 GeV |2j. The resonance /o(1500) shows a flavor-blind decay pattern and is narrow, thus representing 
an optimal candidate to be (predominantly) a scalar glueball. It has been investigated in a large variety of works, e.g. 
Refs. 0,0] and refs. therein, in which mixing scenarios involving the scalar resonances /o(1370), /o(1500), and /o(1710) 
are considered. The second lightest lattice-predicted glueball state has tensor quantum numbers (J PC = 2 ++ ) and a 
$_i | mass of about 2.2 GeV; a good candidate could be the very narrow resonance /j(2200) @, if the total spin of the 
latter will be experimentally confirmed to be J = 2. 

The third least massive glueball predicted by lattice QCD has pseudoscalar quantum numbers (J PC = '") and a 
mass of about 2.6 GeV. Quite remarkably, most theoretical works investigating the pseudoscalar glueball did not take 
into account this prediction of Yang-Mills lattice studies, but concentrated their search around 1.5 GeV in connection 
with the isoscalar-pseudoscalar resonances 7/(1295), 77(1405), and 7/(1475). A candidate for a predominantly light 
pseudoscalar glueball is the middle-lying state 77(1405) due to the fact that it is largely produced in (gluon-rich) J/tp 
radiative decays and is missing in 77 reactions [9J. In this framework the resonances 7/(1295) and 7/(1475) represent 
radial excitations of the resonances 77 and 77'. Indeed, in relation to 77 and 77', a lot of work has been done in determining 
the gluonic amount of their wave functions. The KLOE Collaboration found that the pseudoscalar glueball fraction 
in the mixing of the pseudoscalar- isoscalar states 77 and 77' can be large (~ 14%) 10]. However, the theoretical work 
of Ref. [ll| found that the glueball amount in 77 and 77' is compatible with zero. 

In this work we study the decay properties of a pseudoscalar glueball state whose mass lies, in agreement with 
lattice QCD, between 2 and 3 GeV. Following Ref. [L2j we write down an effective chiral Lagrangian which couples 
the pseudoscalar glueball field (denoted as G) to scalar and pseudoscalar mesons. We can thus evaluate the widths for 
the decays G — > PPP and G — > PS, where P and S stand for pseudoscalar and scalar quark-antiquark states. The 
pseudoscalar state P refers to the well-known light pseudoscalars {it, K,rj,T]'}, while the scalar state S refers to the 
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quark- ant iquark nonet of scalars above 1 GeV: {ao(1450), Kq (1430), /o(1370), /o(1710)}. The reason for the latter 
assignment is a growing consensus that the chiral partners of the pseudoscalar states should not be identified with 
the resonances below 1 GeV, see Refs. [IMEl for results within the so-called extended linear sigma model and also 
other theoretical works in Ref. [|[ (and refs. therein). 

The chiral Lagrangian that we construct contains one unknown coupling constant which cannot be determined 
without experimental data. However, the branching ratios can be unambiguously calculated and may represent a 
useful guideline for experimental search of the pseudoscalar glueball in the energy region between 2 to 3 GeV. In 
this respect, the planned PANDA experiment at the FAIR facility [15j will be capable to scan the mass region above 
2.5 GeV. The experiment is based on proton-antiproton scattering, thus the pseudoscalar glueball G can be directly 
produced as an intermediate state. We shall therefore present our results for the branching ratios for a putative 
pseudoscalar glueball with a mass of 2.6 GeV. 

On the other hand, it is also possible that the pseudoscalar glueball G has a mass that is a bit lower than the lattice 
QCD prediction and that it has been already observed in the BESIII experiment where pseudoscalar resonances have 
been investigated in J/ip decays [HI]. In particular, the resonance V(2370) which has been clearly observed in the 
7r + 7r~?/ channel represents a good candidate, because it is quite narrow (~ 80 MeV) and its mass lies just below the 
lattice QCD prediction. For this reason we repeat our calculation for a pseudoscalar glueball mass of 2.37 GeV, and 
thus make predictions for the resonance V(2370), which can be tested in the near future. 

This paper is organized as follows. In Sec. II we present the effective Lagrangian coupling the pseudoscalar glueball 
to scalar and pseudoscalar quark-antiquark degrees of freedom, and we calculate the branching ratios for the decays 
into PPP and SP. Finally, in Sec. Ill we present our conclusions and an outlook. 



II. THE EFFECTIVE LAGRANGIAN 

Following Ref. [l2| we introduce a chiral Lagrangian which couples the pseudoscalar glueball G = \gg) with quantum 
numbers J PC = *" to scalar and pseudoscalar mesons 

= ic d<i ,G (det* - det$ f ) , (1) 

where Cq^ is a coupling constant, 

$ = (S a + iP a )t a (2) 

represents the multiplet of scalar and pseudoscalar quark-antiquark states, and t a are the generators of the group 
U(Nf). In this work we consider the case Nj — 3 and the explicit representation of the scalar and pseudoscalar 
mesons reads LLJ, LL7fl : 



1 

V2 



a^+in- -*°) K o +iK o j . (3) 

V K s + iK~ K° s + iK° 




Under Ul(3) x Ur(3) chiral transformations the multiplet $ transforms as $ — > Ul$U r where Ul and Ur are U(3) 
matrices. The determinant of $ is invariant under SU(3)l x SU(3)r, but not under U(1)a- On the other hand, the 
pseudoscalar glueball field G is invariant under U(3)l x J7(3).r transformations. Under parity, <& — > $t and G —> — G, 
thus the effective Lagrangian of Eq. (JTJ is invariant under SU(3)l x SU(3)r and under parity. Notice that Eq. 
(fTJ) is not invariant under C/a(1), in agreement with the so-called axial anomaly in the isoscalar-pseudoscalar sector. 
The rest of the mesonic Lagrangian which describes the interactions of $ and also includes (axial-)vector degrees of 
freedom is presented in Sec. IA II of the Appendix. For more details, see Refs. (l7l Il8j. 

The assignment of the quark-antiquark fields in this paper is as follows: (i) In the pseudoscalar sector the fields tt 
and K represent the pions or the kaons, respectively 0]. The bare fields t]n = \uu + dd) /V2 and 775 = \ss) arc the 
non-strange and strange contributions of the physical states rj and 77' [|[ : 

V = Vn cos <p + rjs sin ip, i]' = sin ip + i] S cos (p, (4) 

where <p ~ —44.6° is the mixing angle [l7j- Using other values for the mixing angle, e.g. ip = —36° [l9| or ip — —41.4°, 
as determined by the KLOE Collaboration [10| . affects the presented results only marginally, (ii) In the scalar 
sector we assign the field So to the physical isotriplet state ao(1450) and the scalar kaon fields K$ to the resonance 
i4Tg(1430). Finally, the non-strange and strange bare fields on = \uu + dd) and as = \ss) are assigned to the 
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physical isoscalar resonances /o(1370) and /o(1710). The mixing of the bare fields a^ and as is small [3 53; (3 ( m 
agreement with large- N c arguments) and is neglected in this work, but should be taken into account in more refined 
calculations which include, in addition to scalar quark- antiquark states, also the scalar glueball. 

In order to evaluate the decays of the pseudoscalar glueball G we have to take into account that the spontaneous 
breaking of chiral symmetry takes place, which implies the need of shifting the scalar-isoscalar fields by their vacuum 
expectation values 4>N and 0s, 

a N -> a N + cj) N and a s -t a s + <j>s ■ (5) 

In addition, when (axial-)vector mesons are present in the Lagrangian, one also has to 'shift' the axial-vector fields 
and to define the wave-function renormalization constants of the pseudoscalar fields: 

7? -> Z^tt , K l -> Z K K\ rjj -> Z ni r\j , (6) 

whereas i — 1,2,3,4 runs over the four kaonic fields and j = N,S. The numerical values of the renormalization 
constants are Z„ = 1.709, Zk = 1.604, Zk s = 1.001, Z VN = Z w , Z ns = 1.539 [l7| . Moreover, the condensates and 
4>s read 

(j) N = z„f v = 0.158 GeV, <j) S = 2Zk ^ k _^ ^ n = 0.138 GeV , (7) 

v 2 

where the standard values / w = 0.0922 GeV and Jk — 0.110 GeV have been used Q. Once the operations in Eqs. 
(0 and ([5]) have been performed, the Lagrangian in Eq. (TTJ) contains the relevant tree-level vertices for the decay 
processes of G, see Appendix (Sec. IA2| ). 

The branching ratios of G for the decays into three pseudoscalar mesons are reported in Table I for both choices 
of the pseudoscalar masses, 2.6 and 2.37 GeV (relevant for PANDA and BESIII experiments, respectively). The 
branching ratios are presented relative to the total decay width of the pseudoscalar glueball T*?*. (For details of the 
calculation of the three-body decay we refer to Sec. IA 31 of the Appendix.) 



Quantity 


Case (i): M 6 = 2.6 GeV 


Case (ii): M 6 = 2.37 GeV 


1 G-+KKv' G 


0.049 


0.042 


1 G^-KKv' 1 G 


0.019 


0.011 


r - /r*- ot 

1 G—>vvv' G 


0.016 


0.013 


G^vvv'/ 1 G 


0.0017 


0.00080 


p /pTOC 

G->vv'v'' L G 


0.00013 





/T^tOt 

1 G^KKttI l g 


0.46 


0.46 


-p _ /-ptot 
1 G-iriTTTl G 


0.16 


0.16 


P ^ /ptOt 

1 G^v'nn/ 1 G 


0.094 


0.088 



TABLE I: Branching ratios for the decay of the pseudoscalar glueball G into three pseudoscalar mesons. 



As a next step we turn to the decay process G — > PS. The results, for both choices of Mq, are reported in Table 
II. Future works in this direction should also take the scalar glueball [and its mixing with quark-antiquark states 
into account, thus allowing to evaluate the decays into the three scalar-isoscalar resonances /o(1370), /o(1500), and 
/o(1710). 



Quantity 


Case (i): M 6 = 2.6 GeV 


Case (ii): M e = 2.37 GeV 


p /p^ot 
A G^KKs' G 


0.059 


0.069 


p /ptOt 

1 G->-an7r/ A G 


0.083 


0.10 


p /ptOt 

1 G-» 7 <T N / i G 


0.028 


0.033 


p /ptOt 

1 G-lws/ 1 G 


0.012 


0.0092 


p /ptot 
1 G->J7'o- N / 1 G 


0.019 


0.013 



TABLE II: Branching ratios for the decay of the pseudoscalar glueball G into a scalar and a pseudoscalar meson. 
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FIG. 1: Solid (blue) line: Total decay width of the pseudoscalar glueball with the bare mass Mq= 2.6 GeV as function of the 
coupling Cq^. Dashed (red) line: Same curve for Mq= 2.37 GeV. 



In Fig. [T] we show the behavior of the total decay width T 1 ? 1 — F G ^ ppp + F G ^ PS as function of the coupling 
constant c G<s , for both choices of the pseudoscalar glueball mass. (We assume here that other decay channels, such 
as decays into vector mesons or baryons are negligible). In the case of M G = 2.6 GeV, one expects from large- iV c 
considerations that the total decay width T^* < 100 MeV. In fact, as discussed in the Introduction, the scalar glueball 
candidate /o(1500) is roughly 100 MeV broad and the tensor candidate /j(2220) is even narrower. In the present 
work, the condition T*?* < 100 MeV implies that c G ^ < 5. Moreover, in the case of M G = 2.37 GeV in which the 

identification G = X(2370) has been made, we can indeed use the experimental knowledge on the full decay width 
[ r x(2370) = 83 ± 17 MeV [l6|] to determine the coupling constant to be c G ^ = 4.44 ± 0.45. 
Some comments are in order: 

(i) The results depend only slightly on the glueball mass, thus the two columns of Table I and II are similar. It 
turns out that the channel KKtt is the dominant one (almost 50%). Also the r/inr and ri'irir channels are sizable. On 
the contrary, the two-body decays are subdominant and reach only 20% of the full mesonic decay width. 

(ii) Once the shifts of the scalar fields have been performed, there are also bilinear mixing terms of the form GrjN 
and Grjs which lead to a non-diagonal mass matrix. In principle, one should take these terms into account, in addition 
to the already mentioned t/nVS mixing, and solve a three-state mixing problem in order to determine the masses of 
the pseudoscalar particles. This will also affect the calculation of the decay widths. However, due to the large mass 
difference of the bare glueball fields G to the other quark-antiquark pseudoscalar fields, the mixing of G turns out to 
be very small in the present work, and can be safely neglected. For instance, it turns out that the mass of the mixed 
state which is predominantly glueball is (at most) just 0.002 GeV larger than the bare mass M G = 2.6 GeV. 

(iii) If a standard linear sigma model without (axial-)vector mesons is studied, the replacements Z„ = Zj^ = Z VN = 
Z„ s — 1 need to be performed. Most of the results of the branching ratios for the three-body decay are qualitatively, 
but not quantitatively, similar to the values of Table I (variations of about 20%). However, the branching ratios for 
the two-body decay change sizably w.r.t. the results of Table 2. This fact shows once more that the inclusion of 
(axial-)vector degrees of freedom has sizable effects also concerning the decays of the pseudoscalar glueball. 

(iv) In principle, the three-body final states for the decays shown in Table I can also be reached through a sequential 
decay from the two-body final states shown in Table II, where the scalar particle S further decays into PP, for 
instance, Kq(1A30) — > Kir. There are then two possible decay amplitudes, one from the direct three-body decay and 
one from the sequential decay, which have to be added coherently before taking the modulus square to obtain the 
total three-body decay width. Summing the results shown in Table I and II gives a first estimate (which neglects 
interference terms) for the magnitude of the total three-body decay width. We have verified that the correction 
from the interference term to this total three-body decay width in a given channel is at most of the order of 10% 
for M G = 2.6 GeV and 15% for M G = 2.37 GeV. For a full understanding of the contribution of the various decay 
amplitudes to the final three-body state, one needs to perform a detailed study of the Dalitz plot for the three-body 
decay. 
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III. CONCLUSIONS AND OUTLOOK 



In this work we have presented a chirally invariant effective Lagrangian describing the interaction of the pseudoscalar 
glueball with scalar and pseudoscalar mesons for the three-flavor case Nf = 3. We have studied the decays of the 
pseudoscalar glueball into three pseudoscalar and into a scalar and pseudoscalar quark-antiquark fields. 

The branching ratios are parameter-free once the mass of the glueball has been fixed. We have considered two 
possibilities: (i) in agreement with lattice QCD we have chosen Mq = 2.6 GeV. The existence and the d ecay 
properties of such a hypothetical pseudoscalar resonance can be tested in the upcoming PANDA experiment [15]. 
(ii) We assumed that the resonance A (2370), measured in the experiment BESIII, is (predominantly) a pseudoscalar 
glueball state, and thus we have also used a mass of 2.37 GeV [16j . The results for both possibilities have been 
summarized in Tables I and II: we predict that KKn is the dominant decay channel, followed by (almost equally 
large) rjTrir and r/7T7T decay channels. In the case of BESIII, a measurement of the branching ratio for other decay 
channels than the measured r/'mr one could ascertain if A(2370) is (predominantly) a pseudoscalar glueball. In the 
case of PANDA, our results may represent a useful guideline for the search of the pseudoscalar glueball. 

Future studies should consider possible mixing of the pseudoscalar glueball with charmonia states and include the 
light scalar glueball in the scalar sector. Moreover, the mechanism of the glueball production via proton-antiproton 
fusion can also be investigated. 
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Appendix A: Details of the calculation 



1. The full mesonic Lagrangian 



The chirally invariant U(Nf)L x U(Nj)n Lagrangian for the low-lying mesonic states with (pseudo)scalar and 
(axial-)vector quantum numbers has the form 

C mes = ■&[(£>„$)+(£>„*)] - m2Tr($t$) _ Aifft(*t$)]3 _ A 2 Tr($t$)2 



and 



-Tr[(Ln 2 + (Rn 2 } + Tr[(^i + A)(L* + + Tr[ff($ + $t)] 
ci(det$ - dct$t) 2 + i^{Tr(L llv [L' 1 ,L v ]) + Tr(R^[R^ R"})} 
- ^Tr($t$ )Tr ( L 2 + + ^TrOZ^I 2 + \<S>R,\ 2 } 
■2h 3 Tr(L li &Ri i &). 



(Al) 



where 



L u = 



1 

V2 



V2 



V2 



«>N-P' 



•J2 _V% 



K*^+ + iff + \ 



— K*^° + K 



no 



7& J 



R,i = 



1 

V2 



V2 



K*po 



J IN °1 

_v / 2 



K 



Jli 



for details see Refs. 
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In the present context we are interested in the wave- function renormalization constants Z^ introduced in Eq. ([6]). 
Their explicit expressions read [l?], EH : 



Z« = Z m = ai , (A2) 

'm 2 ai - gl4> 2 N 



Z K = 2mK > , (A3) 



Zk s = ^ , (A4) 

yj4m 2 K * - gf(<PN - V2<f>s) 2 

Z vs = ; mhs . (A5) 

m /i S - 2 Si^is 



2. Explicit form of the Lagrangian in Eq. ([TJ) 

After performing the field transformations in Eqs. and ©, the effective Lagrangian (JXJ) takes the form: 
£*f = ^G(V2Z Ks Z K a K s K° + V2Z K ZKsa K K° s - 2Z Ks Z K a+ K S R- 

- 2Z Ks Z K a+K s K° - 2Z Ks Z K a^s K+ ~ V2Z Ks Z K a" Kg K+ - V2Z 2 K Z m K a K° m 
+ V2Z] (s Z VN K^K° s r] N - y/2Z%Z m K~ K + rj N + Z vs a^rj S + 2Z ns a Q a^rj S 

+ Zl N Z Vs r) 2 N iis - V2Z%Z n K°K°ir° + ^f2Z 2 Ks Z,K%^ s ^ + V2Z 2 K Z^R- K+tt° 

- Z vs ZlT)sTr° 2 + 2Z] ( ZJ( K + it- + 2Z 2 k Z v K q K->k + - 2Z 2 Ks Z 7r K Q s K s n + 

- 2Z vs Z 2 r]sTr-TT+ - 2Z Ks Z K a^K+K° + V2Z 2 Ks Z m K+ K sVn - V^Z^Z^ K^t: 

- 2Z 2 Ks Z 7T K+K° S 7r- - V2Z Ks Z K a° K+K- + V2Z K Z Ks K- K+(j> N + V2Z K Z Ks R- K+a N 
+ V2Z Ks Z K K° s K l Vat + V2Z Ks Z K K s K°a N + V2Z K Z K3 K K° S (f> N + V2Z K Z Ks K°Tf s a N 
+ V2Z Ks Z K KgK + (f) N + V2Z Ks Z K Kg K + a N - Z ns T] S (f> 2 N - Z ns ris<7 2 N - 2Z r]s r\ s 4>NON 

+ 2Z 7r a,Qn (f) S + 2Z 1T aQTT°<rs + 2ZnO>o ir~<t>s + 2Z ff a^"7r~ <t s + 2Z 7T a^ ) tt + (f> s + 2Z 7T a^ ) ir + <r s 

- 2^ j?n ?7jv'/ , jv0s - ^Z^ N r\ N <i> N a s - 2Z riN rj N a N (/)s - 2Z rw r\ N a N a s ) . (A6) 

The latter expression is used to determine the coupling of the field G to scalar and pseudoscalar mesons. 

3. Tree-body decay 



For completeness we report the explicit expression for the three-body decay width for the process G — > P\PiP; 
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32(27r) 3 M| 7(„ ll+m2 )2 J{m 2a )n 

where 



(m 23 ) mi n = (i? 2 * + ^) 2 - ( ^jEf-mj + jEf-m 2 ) , (A7) 



(m 23 ) max = (E* 2 + El) 2 - (jEf-m 2 - ' Ef - mU 



2 



(A8) 
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and 

^ ~ 2^ ' E * ~ 2^ ■ (A9) 

The quantities mi, m2, TO3 refer to the masses of the three pseudoscalar states Pi, P2, and P3, -Mg^p p p i s the 
corresponding tree-level decay amplitude, and s g->p 1 p 2 p 3 is a symmetrization factor (it equals 1 if all Pi, P2, and P3 
are different, it equals 2 for two identical particles in the final state, and it equals 6 for three identical particles in the 
final state). 

For instance, in the case G — > K~K + tt° one has: | — i-M-Q^K- k+tt° I 2 = \ C %^K^' 1711 = 1712 = mK = 0-494 
GeV, m 3 = m w o = 0.135 GeV, and M G = 2.6 GeV. Then: 

r G ^ K -K+*o = 0.00041 c| $ [GeV] . (A10) 
The full decay width into the channel KKtt results from the sum 

^G^KKtt = ^G^K- K+ir° + ^G^K°K °7r° + ^G^K°K+tt- + ^G^K°K-7r+ = ^G^K~K+w° 1 (AH) 

The other decay channels can be calculated in a similar way. 
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